Astron. Nachr. / AN 334, No. 1, 1 -|7](2013) / DOI This. is/not. aDOI 



The stratified evolution of a cool star 
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A low mass star usually experiences stratification and abundance anomalies during its evolution. A 0.95 Mq star with a 
metallicity Z — 0.004 is followed from the main-sequence to the Horizontal Branch (HB). On the main-sequence the 
larger effects of stratification may come from accretion as was suggested in relation to metallicity and planet formation. As 
it evolves through the giant branch, stratification appears around the hydrogen burning shell. It may create hydrodynamic 
instabilities and be related to abundance anomalies on the giant branch. After the He flash the star evolves to the HB. If it 
loses enough mass, it ends up a hot HB star (or in the field an sdB star) with effective temperatures larger than 1 1000 K. 
All sdB stars are observed to have an approximately solar iron abundance whatever their original metallicity, implying 
overabundances by factors of up to 100. So should the 0.95 solar mass star. How its internal hydrodynamic properties on 
the main sequence may influence its fate on the HB is currently uncertain. 



1 — '■ 1 Astrophysical context 



In most stellar mass intervals, stars are significantly strat- 
ified by particle transport processes during part of their 
evolution. This is in addition to stratification due to nu- 
clear burning. For instance, as a 0.95 Mq star of metallic- 
ity Z — 0.004 remains cooler than 6500 K on the main- 
sequence, the larger effects of potential stratification come 
possibly from accretion but are limited by thermohaline 
mixing. This could be linked to planet accretion and the de- 
struction of Lithium. 

As it ascends the giant branch, effects are very small at 
the surface but might be larger close to the H burning shell. 
The largest stratification occurs on the horizontal branch 
(HB) but is strongly dependent on the effective temperature 
that the star ends up with. The spread of HB masses that 
leads to the effective temperature spread on the HB is po- 
tentially linked to the low mass star structure on the main- 
sequence. 

We here review the stratification during the evolution 
of such a star using mainly results obtained with an evolu- 
tionary code that takes into account transport by atomic dif- 
fusion outside of convection zones. Calculations were car- 
ried out with all usual equations of standard stellar evolution 
with in addition a set of 56 coupled differential equations to 
take into full acount the effects of the atomic diffuson of the 
28 species included. Radiative acceleration and Rosseland 
averaged opacity are continuously recalculated during evo- 
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lution so that the calculations are continuously consistent 
with all abundance variations as described in Richer et al. 
d!998t and Turcotte et al. J 1998b where the only parameter 
of the calculations, the mixing length, is fixed by proper- 
ties of the Sun. Parameters are introduced only to describe 
macroscopic physical processes competing with atomic dif- 
fusion. 

The structure of this review follows the major evolution- 
ary stages of a low mass star; we are not concerned with 
stratification originating directly from nuclear reactions. It 
is stratification caused by transport processes that is looked 
into here. 

2 Main sequence evolution of a cool star 

A 0.95 M Q star of original metallicity Z = 0.004 was fol- 
lowed from the pre— main— sequence to the He flash. As 
may be seen from Fig. 1 of Michaud, Richer & Richard 
( 120101 1 the effects of diffusion on the HR diagram become 
apparent mainly around turnoff as well as close to the hook, 
where on the giant branch the H burning shell crosses the 
point of deepest inward expansion of the surface convection 
zone. These effects are small. 

The variations of the surface abundances of He and 
Fe are shown in Fig[T] They were calculated for all 28 
species but the reduction factors are similar for most met- 
als. The depth of the convection zone causes the radiative 
accelerations to play a small role which may be seen from 
Fig.[2]where the interior concentrations for all 28 calculated 
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Fig. 1 Surface abundances of He and Fe during the evolu- 
tion of a 0.95 Mq star from the ZAMS to the He flash. The 
settling is progressive during evolution and the underabun- 
dance is by the largest factor (-0.15 dex for He) at turnoff. 
The He abundance returns to its original value as the He 
flash is approached but the Fe abundance remains -0.01 dex 
below as some Fe is concentrated in the stellar center. 



species are shown at an age of 8.6 Gyr. The interior con- 
centration variations caused by diffusion processes are by 
around 0. 1 dex (note that the concentration scale on the right 
inset of the figure ranges from a factor of 0.85 to 0.95). For 
light elements (from H to O) nuclear raections are the main 
cause of concentration changes and are concentrated in the 
inner 50% by radius. In the outer 50% of the radius the 
changes are due to atomic diffusion which is here domi- 
nated by gravitational settling. However one sees the effect 
of radiative accelerations below the convection zone from K 
to Mn. 



2.1 Metallicity around 13 Gyr 

The Pop II stars with the lowest metallicity (see for instance 
Beers & Christlieb 2005 ) are thought to have formed some 
13 Gyr ago, a few hundred Myr after the Big Bang. They 
may be used as messengers from the earliest times. But to 
what extent are we certain that their surface composition 
we see today is that with which they formed. The 0.95 Mq 
star we follow will not see its surface abundances modi- 
fied strongly by atomic diffusion. But this is not generally 
the case for halo fields stars as may be seen in Fig. [3] As 
one considers stars with the lowest metallicities, some of 
the surviving ones may have effective temperatures larger 
than 5800 K. If there is no process outside of convection 
zones competing wih atomic diffusion, surface overabun- 
dances of iron by more than a factor of ten are possible so 
that the apparent metallicity would be larger than the one 
the star formed with. This could be detected since iron to 
carbon or to oxygen ratios would appear much larger than 
solar since the carbon and oxygen abundances are not ex- 
pected to be larger than the original one (see Fig. 9, 10 and 
11 of Richard, Michaud & Richer 120021 . Most but not all 
observed very low metallicity stars have an effective tem- 
perature smaller than 5800 K. 




Fig. 2 The interior concetrations of the 28 calculated 
species are shown intensity coded at an age of 8.6 Gyr. The 
radial coordinate is the radius and its scale is linear, but the 
logarithmic value of the mass coordinate above a number 
of points, log AM/M t , is shown on the left of the horizon- 
tal black line. The concentration scale is given in the right 
insert. Small circles mark the surface convection zone. The 
effect of radiative accelerations is evident from K to Mn 
slightly below the convection zone. It is small as will be the 
case throughout the main— sequence. 

2.2 Accretion of planets 

There is now ample evidence for the existence of planets 
around most stars. Through disk migration and tidal inter- 
actions, planets move inwards leading to accretion (e. g. 
Jackson, Barnes & Greenberg [20091 ) ■ Many of the observed 
innermost exoplanets are expected to be accreted within a 
few Gyrs. Since planets have larger metallicity than the stel- 
lar cloud from which they formed, they should increase the 
metallicity of the accreting star. On the other hand stars 
that have been observed to have planets have been found to 
have larger metallicity than those without planets (Fischer 
& Valenti 2005). Is this overmetallicity of stars with planets 
due to the accretion of planets or is an overmetallicity of the 
original nebula from which the star and planet formed the 
cause of planet formation? In other words, is the overmetal- 
licity primordial or due to planet accretion? The answer to 
this question depends on the mass in which the planet mixes 
after merging with the star. 

Vauclair (2004) suggested that the larger metallicity of 
the matter accreted from the planet implied that, due to 
thermohaline mixing, the mass into which the planet mixed 
was generally much larger than the mass of the convection 
zone. Consider a planet with a 1.0 Jupiter mass accreting 
to a 1.4M© star. Garaud (120111) found that (see her Fig. 2) 
within 10 4 yr, the depth of the convection zone adjusts so 
that the Ledoux criterion is satisfied at its bottom; over- 
metallicity has been reduced by 20%. Thermohaline con- 
vection further, but more slowly, extends the mixing. After 
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Fig. 3 The colored line segments represent the iron abun- 
dance between 12.6 and 13.2 Gyr in stars of different mass 
and different original metallicity. Stars which are before 
turnoff appear as a point, since slow evolution causes the 
segment to collapse to a point. However those which are 
in more rapid evolutionary stages (around turnoff or on the 
subgiant branch) appear as line segments. Note that subgiant 
branch stars have slightly larger abundances than main- 
sequence stars of the same T e ff. The initial Fe abundance 
is indicated by a small line on the left. 

10 7 yr, overmetallicity has been reduced by a factor of ten 
so that the effect of the accretion has virtually disappeared 
from the surface. Consequently, the increased metallicity 
observed on planet bearing stars has to be of primordial ori- 
gin. It cannot have been caused by planet accretion since the 
effect of the accretion lasts for only a small fraction of the 
star's evolution time. Planet accretion cannot be an impor- 
tant cause of star stratification (see also Proffitt & Michaud 
119891 ). These results depend, however, on an evaluation of 
the efficiency of thermohaline convection. 

2.3 Thermohaline convection 

Thermohaline convection is for instance observed in 
oceanic waters when there is salt water above cooler fresh 
water. There it is often called salt fingering. In stars the in- 
creased metallicity plays the role of the increased salinity. 
In the absence of a proper calculation of the transport ef- 
ficiency of thermohaline convection in stars, one has often 
used a mixing length evaluation. It is however very uncer- 
tain since it depends rather sensitively on the assumed elon- 
gation of the cells. Different evaluations differed by a factor 
of ~ 100. More recently detailed simulations first in 2-D 
(Denissenkov 2010 ) and then in 3-D (Traxler, Garaud & 



Stellmach l201 II) have improved the situation. The results of 
the two simulations are in agreement. 

In their 3-D simulation, Traxler, Garaud & Stellmach 
( 1201 11 1 determine the turbulent diffusivity in a box where 
they impose a stably stratified thermal field but an unstably 
startified compositional field. They carried out simulations 
at resolutions of 96 x 96 x 192 to 384 x 384 x 384 in order to 
approach as much as possible the parameters appropriate for 
the stellar case. While they were still far from it, they claim 
to have arrived at an asymptotic value for the turbulent dif- 
fusivity that may be used as a universal law. As shown in 
their Fig. 3, they agree with the evaluation of Kippenhahn, 
Ruschenplatt and Thomas ( 119801 ) which is the lower of the 
evaluations previously considered. While their results seem 
to converge to an asymptote as the ratio of thermal to com- 
positional diffusivity and the Prandtl number approach stel- 
lar values, they are still far from them and higher resolution 
simulations would still appear needed to firmly establish the 
results. Their results were used in the discussion of § !2.2l and 
have some impact also on the transport above the H burning 
shell of giants described in the next section as well as for 
the Li abundance in stars (Theado & Vauclau T2012l ). 

3 Giant Branch 

As the 0.95 M Q , Z = 0.004 star goes up the giant branch, 
its hydrogen burning shell approaches the depth where the 
surface convection zone had its deepest extension (at the 
first dredge-up). Concentrations were homogenized from 
the surface down to that depth but there appeared a disconti- 
nuity in many concentrations at that depth. As the hydrogen 
burning shell crosses that discontinuity, the star describes 
what is often called the hook in the HR diagram. As the 
hook is approached there occurs, slightly ahead of the H 
burning shell, an H abundance increase caused by 3 Hc burn- 
ing (the dotted gray line on the lower panel of Fig.|U). This 
H abundance translates into a /i gradient inversion, shown 
as the dotted gray line in the upper panel of Fig.|4l which 
has been suggested to cause mixing, by thermohaline con- 
vection, between the burning shell and the surface (Eggle- 
ton, Dearborn & Lattanzio l2006l Charbonnel & Zahn l2007t . 
This could explain abundance variations seen on the giant 
branch. 

However atomic diffusion starts to affect the homoge- 
nized concentrations as soon as the convection zone starts 
receding. Does atomic diffusion have a significant effect on 
the fi gradient inversion that drives the thermohaline con- 
vection? 

On the lower panel of Fig.|4]are shown (from data used 
by Michaud, Richer & Richard (2010) where more details 
may be found; see in particular their Fig. 8) the X(H) in- 
terior profiles immediately after the convection zone starts 
receding (black lines) and just before the H burning shell 
crosses the region of the compositional break (gray lines). 
Solid lines represent calculations with diffusion and dotted 
lines those without. The H abundance change is ~ 3 times 
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Fig. 4 Upper panel Interior profile of p as a function of 
the interior mass, M r , immediately after the first dredge- 
up (black lines) and (gray lines) as the hydrogen burning 
shell crosses the region of the composition break left by 
the first dredge-up (more precisely ~ 1.5 Myr before). The 
solid lines represent the calculations with atomic diffusion 
and the dotted lines those without. The solid lines refer to 
the scale on the left and the dotted lines to the scale on the 
right (the two scales were adjusted so that the two horizontal 
black lines are superposed). The two scales cover the same 
A(l//i) interval. Lower panel Interior profiles of the corre- 
sponding X(H) values. 

larger in the model with diffusion than in that without diffu- 
sion while the effect on p is by a factor of ~ 6. It would then 
appear that the gravitational settling of He leading to the H 
abundance increase has a larger effect than 3 He burning. 

The turbulence induced by thermohaline convection 
tends to counteract gravitational settling. The concentration 
gradients of metals and of He (or H) implied by the val- 
ues of X(H) on Fig.|4]are very small, yet it is only through 
those gradients that turbulence (Z>r) has an effect, whereas 
atomic diffusion also acts through the much larger g driv- 
ing terms. Consequently, Di can have an effect only if it is 
much larger than the atomic diffusion coefficient, Di p . But 
how much larger? To have a significant effect on the diffu- 
sion velocity, Dj must lead to a contribution similar to the 
driving terms of atomic diffusion. Turbulence then has an 
effect if 



D 1 



di 



D 



kT 



-f) 



(1) 



Eq. (Q]) was evaluated roughly using drift velocities of He 
and the value of X(He) from Michaud et al ( |20T0] l. The D T 
so obtained from Eq. (Q]) is shown on Fig. 9 of that paper. It 



varies between 10 6 , immediately above the burning shell, 
and 10 7 cm 2 /s immediately below the surface convection 
zone. According to Sect. 12. 31 this upper limit is two orders 
of magnitude larger than the best evaluations of the thermo- 
haline convection turbulent diffusion coefficient. The gra- 
dients that atomic diffusion leads to should then easily be 
maintained in the presence of this turbulence. 

Whether the p inversion has any importance for surface 
abundances also depends on the mixing coefficient that ther- 
mohaline convection leads to. The lower value of the mixing 
coefficient discussed in Sect. l2.3l is below what is required 
to cause sufficient mixing with the surface according to Ga- 
raud ( 201 1 ). That conclusion is however based only on the 
contribution of 3 He burning to the p gradient inversion since 
atomic diffusion was not included in the calculations she 
refers to. Would the larger p gradient values that atomic dif- 
fusion leads to make a difference? A calculation including 
atomic diffusion and thermohaline convection with the co- 
efficient described in Sect. l2.3l is needed to settle this ques- 
tion. 



4 Horizontal branch and sdB stars 

As a 0.95 M star leaves the giant branch, it proceeds to the 
horizontal branch (HB). Some such stars loose more mass 
than others through a process which is not currently well 
understood. Depending on the mass it looses the star ends 
up with a different T e ff on the HB. In Pop II stars, some HBs 
of globular clusters are populated to relatively high T e g and 
the field equivalent, the sdB stars, extend into the sdO stars. 

In globular clusters, HB stars should have the same com- 
position as on the giant branch of the cluster except perhaps 
for species involved in the CNO cycle and affected by pro- 
ton burning. The abundance of most metals, and in partic- 
ular of iron, should be the same as on giants. However it 
had been predicted that atomic diffusion should cause abun- 
dance anomalies on the HB (Michaud, Vauclair & Vauclair 
|19831 >. This prediction has now been strikingly confirmed in 
all globular clusters with a HB extending to T e ff > 1 1000 K. 
As shown on Fig.[5]for the clusters M15, M68 and M92, HB 
stars hotter than 11000K have approximately a solar iron 
abundance while the cooler ones have the same Fe abun- 
dance as giant stars of the cluster (~ 1/200 solar). Simi- 
larly in clusters with different metallicities (e.g. M3, M13, 
NGC288, NGC6752, NGC1904, NGC2808), HB stars with 
T e ff > 11000K have about a solar iron abundance while 
the lower T e s stars have the same iron abundance as giants 
(Behr|2003 Moehler et al [20001 Fabbian et al 120051 Pace 
et al |2006l >. The black points on Fig. [5] show a few excep- 
tions: a few more rapidly rotating stars around 1 1000 K. The 
black line segments show the calculated surface abundances 
of iron in 8 HB models covering the T e ff range. The dotted 
part covers the first 10 Myr of HB evolution and the solid 
part, the following 20 Myr. The observed Fe abundances 
are as expected for most stars with T e ff > 11000 K. For the 
cooler stars and for the three represented by black points 
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Fig. 5 Upper panel Observed Fe abundances in HB stars 
of M15, M68 and M92. The black line segments represent 
the iron abundances during the HB. Dotted part, first 10 Myr 
and solid part from 10 to 30 Myr. Lower panel Correspond- 
ing rotational velocities. Gray points are data from Behr 
(12003b while the three black ones identify rapidly rotating 
ones. 

around 1 1000 K, rotation (lower panel of Fig-Eb appears to 
compete succesfully with the diffusion driven by radiative 
forces. It has been suggested that, given the observed v sin i 
on the HB, meridional circulation is efficient in reducing 
abundance anomalies in stars cooler than 10000K but not 
in those hotter than 1 1000K with a buffer zone in between 
(Quievy et al l2009l . 

The calculations used in Fig. [5] are the continuation 
of those described in the preceding sections. In pursuing 
to the HB, approximately the same procedure as used by 
Sweigart( 1987) was followed. After the He flash, this pro- 
cedure involves removing a varying amount of mass in or- 
der to arrive at different T e ffS on the HB. When on the HB 
the model is reconverged and, in the models described here, 
atomic diffusion is allowed to proceed. In these calculations, 
the outer 10~ 7 M Q was assumed to be mixed by some tur- 
bulent process. This is the one adjusted parameter for these 
calculations and it has the same value for all HB and sdB 
stars discussed here. An example of the internal abundance 
distribution this leads to is shown in Fig. [6] (from Michaud, 
Richer & Richard [2011 al l. The effects of atomic diffusion 
are not only superficial but are felt over the outer third of 
the radius of the star. The abundance variations are large, by 
factors of order ten, and one sees clearly the effect of closed 
electronic shells on radiative accelerations in the shifting in- 
wards of the local maxima of concentrations from Ca to Ni. 

On Fig.|7]the calculated Fe surface abundaces are com- 
pared to the observed values in sdB stars. Some 60 evo- 




Fig. 6 The interior abundances of the 23 species included 
for the HB calculations are shown intensity coded at an age 
of 31.1 Myr after the beginning of the HB in a ~ 25000 K 
HB star. The radial coordinate is the radius and its scale is 
linear; the logarithmic value of the mass coordinate above a 
number of points, log AM/iW*, is shown on the left of the 
horizontal black line. The concentration scale is given in the 
right insert. Small circles mark the central convection zone. 
The effect of radiative accelerations is evident from Ne to 
Ni. Note that the concentration scale is different from that 
used for the main— sequence star (see Fig.O . 
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Fig. 7 Observed Fe abundances in sdB stars (circles) and 
upper limits (downward pointing triangles) from Geier et al 
(12010b who observed 139 sdBs with high resolution spec- 
trographs. All stars were verified to be slowly rotating. The 
line segments are models calculated with original Fe abun- 
dances of [Fe/H] = -2.3, -1.3, -0.4 and 0.0. The dotted part 
of the segment is from to 10 Myr after the beginning of 
the HB while the solid segment is from 10 to 32 Myr. The 
red arrow shows the position of a star originally of 0.95 M Q 
with a metallicity Z = 0.004. The masses on the HB are in 
small characters but one may zoom on them in the online 
version. 
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lutionary HB models were calculated with metallicities of 
Z n = 0.0001, 0.001, 0.004 and 0.02. These cover the metal- 
licity range with which one may expect most field sdB stars 
to have formed. During their evolution, even most of those 
starting with a metallicity 200 times smaller than solar end 
up, after 10 Myr, within a factor of ten of solar (the solid part 
of the line segments). This is also the interval were most ob- 
served values are found. The evolutionary calculations were 
stopped after some 32 Myr on the HB for technical reasons. 
The observed sdB stars beyond that time were identified in 
Fig. 2 of Michaud et al (1201 lab by red circles and these are 
also used to identify those objects in the (Fe, T e ff) plane. The 
observed upper limits (the inverted gray triangles) would 
correspond to the first 10 Myr of the HB evolution. 

Similar observations were carried out by the same ob- 
servers for He, C, N, O, Ne, Mg, Al, Si, P, S, Ar, K, Ca, Ti 
and Cr and compared with model calculations with a similar 
level of agreement as for Fe (see Figs. 15 and 16 of Michaud 
et al l2011at . 

5 Conclusion: what do we learn from HB 
stars? 

As seen in the preceding Section, complete evolution- 
ary models taking atomic diffusion and radiative accelera- 
tions into account predict approximately the observed metal 
abundances on hot HB and on sdB stars. Only one param- 
eter was adjusted and the same value of the parameter was 
used for all calculations. This confirms that radiative accel- 
erations are the main cause of the abundance anomalies ob- 
served on these objects. What determines the size of the 
anomalies is the mass between the surface and the layer 
where the separation takes place. This mass could be de- 
termined by some turbulent process and lead to the type 
of solution used here. But when an adjustable parameter 
is involved, in this case the mixed mass, there always re- 
mains the possibility that a different process would have the 
same effect. Indeed, it was found for Pop I stars (Vick et al 
120101 Michaud, Richer & Vick 120 lib! that similar surface 
abundance anomalies are caused at the surface if the pro- 
cess competing with atomic diffusion is either mass loss or 
turbulence so long as the mass loss rate is such that most of 
the chemical separation occurs at the same mass below the 
surface as is mixed in the turbulent model. Can asteroseis- 
mology distinguish between the two possibilities? 

Based on a parameter free equilibrium model, pulsations 
were predicted to occur in some sdB stars (Charpinet et al 
|1996t Fontaine et al 120031 ). They are caused by the Fe ac- 
cumulation where it is the main contributor to opacity. The 
pulsations were observed shortly after being predicted 

The seismic properties of sdB stars were calculated for 
both models with mass loss and models with turbulence as 
the competing process by Hu et al (1201 II) . They conclude 
that the models with turbulence, such as those discussed in 
the preceding sections, ought to be preferred. A more de- 
tailed fit of the pulsating properties of individual stars than 



they attempted would seem to be needed before this conclu- 
sion is definitively accepted. 

Two aspects of hot HB and sdB stars remain puzzling: 
what determines the varying amount of mass lost from the 
giant branch to the HB and why do only the hot stars ro- 
tate slowly. Both of these properties might well be linked 
to properties of the stars while on the lower main sequence 
and possibly to the presence of a magnetic field for instance 
(Charbonneau 2005 ). However, onee the existence of slowly 
rotating sdB and hot HB stars is accepted, their abundance 
anomalies are expected. 
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